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SUMMARY 

Th;s  report  presents  a new  model  for  ray -optical  reflection  from  the  rough,  air/material  interface 
of  a randomly  rough  surface.  We  review  the  current  understanding  of  light  reflection  from  rough  sur- 
faces and  show  how  the  unique  features  of  this  new  model  contribute. 

The  new  model . instead  of  treating  the  rough  > air/material  interface  as  being  composed  of  randomly 
oriented  flat  microareas  (facets),  treats  the  interface  as  an  ensemble  of  randomly  oriented,  randomly 
curved  microareas.  These  two  models  arc  identical  (physically)  and  arc  found  to  give  the  same  results. 
However,  this  new  derivation  leads  to  some  new  and  useful  results.  (DA  new  visualization  of  the  sur- 
face structure  (an  "average”  irregularity)  is  conceived  and  p -wed  valid.  An  average  irregularity  is  an 
optically  smooth  curved  surface  of  revolution  of  a shape  such  that  it  gives  the  'a me  distribution  of 
reflected  light  (when  irudiated  by  a uniform,  well -collimated  beam)  as  that  given  by  the  avtu&!  tough- 
surface  mictostructure  (when  irradiated  by  the  actual,  mmuniform.  well-collimated  beam).  ( ')  It  is 
found  that  the  shape  of  this  surface  of  revolution  mav  he  greatly  restricted  and  still  he  general  enough 
to  represent  any  physically  realistic  microsiructute.  (.  ) Modeling  this  average  irregularity  as  an  ellipsoid 
of  revolution  gives  a surface  structure  function  that  is  more  accurate  and  useful  than  previously  existing 
ones.  (4)  Unlike  the  "facet”  derivations,  this  derivation  lends  itself  to  a normalization  giving  the  abso- 
lute, instead  of  just  the  relative,  reflectance-distribution  function. 

Reflection  bv  this  new  interface  model,  combined  with  some  Lambertian  reflection,  is  tested  ex- 
tensively both  within  and  without  ‘he  nominal  plane  of  incidence  on  a variety  of  commonly  occurring 
rough  surfaces.  This,  we  believe,  is  the  first  time  a iheotv  ui  ray  reflection  from  ihe  interface  is  exten- 
sively tested  outside  the  nominal  plane  of  incidence.  Such  reflection  is  of  particular  interest  to  laser 
target  designation  systems.  Data  for  'M  Black  Velvet  paint,  soiled  olive  drab  paint  cement,  plywood, 
and  grass,  were  used,  and  the  parameters  ui  the  model  were  optimized  to  give  the  best  fit  of  the  model 
to  the  data.  The  comparisons  were  m general  reasonably  good  (rim  deviations  of  .'0  io  .*5  percent),  .md 
the  discrepancies  could  he  mostly  explained  hy  die  existence  ot  gov\s  deviations  of  come  «»i  the  measured 
surlaces  from  the  assumed  surface,  such  as  the  existence  of  two  surface  materials  instead  of  gist  one  am) 
the  existence  of  a significant  ditcclioital  dependence  of  the  surface  structure  Also,  the  normalization  was 
verified  in  a rough  manner. 
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NOMENCLATURE 


Symbol  Definition 

A 

/ Unit  vector  pointing  in  the  direction  from  which  the  light  is  incident. 

r Unit  vector  pointing  in  the  direction  of  reflection. 

n Unit  vector  pointing  in  the  direction  of  the  microsurface  normal. 

0 Zenith  angle  of  incidence,  relative  to  the  macrosurface. 

6,  Elevation  and  azimuth  angles  of  reflection,  relative  to  the  macrosurface. 

a,  a Zenith  :r.i  azimuth  angle  'f  the  micromrfae*  nomtal.  relative  to  the  macrosurface. 

s Zenith  angle  of  incidence  relative  to  the  microsurface 

Angle  between  the  microscopic  and  macroscopic  planes  of  incidence. 

7 Angle  between  the  incident  polarization  plane  and  the  macroscopic  plane  of  incidence. 

n.  k Indices  of  refraction  and  absorption  of  the  surface  material. 

4Us)  Fresnel  reflection  coefficient  of  the  surface  material  (any  subscript  indicates  the  polarization 

state  of  the  incident  light). 

ilA(  A small  area  on  the  microsurface. 

Dftt)  The  surface  structure  function,  or  the  relative  amount  of  micmsurface  area  oriented  in  a 
given  direction. 

c Ratio  of  the  axis  revolved  about  to  the  axis  revolved  for  the  ellipsoid  of  revolution  average 

iurface-irregulaniy 

f>o(ul  Radius  of  eurvatiue  of  the  surface  of  revolution  average  surface -irregularity  in  the  plane 
through  the  rotation  axis. 

p l*,#  Radnu  of  curvature  of  ihe  onfave  of  icsoiuinm  average  surface -irregularity  n.  a plane  per- 
pendicular to  the  rotatHW  axis 

f(l  UR  RDF  ol  bidirectional  reflected -radiance -distribution  fund  too  The  ratio  of  the  radiance 

reflected  (into  a chosen  direction)  lo  the  irradiance  incident  in  a collimated  tseam  (from  a 
chosen  direction). 

f(j  URIDF  or  bidirectional  tcflcctcd-mtcmiis -distribution  function.  The  ratio  of  the  radiant 

intensity  |VV;s»*’|  reflected  (into  a cb>«cn  direction)  to  the  radiant  posset  |\V|  incident 
(from  a chosen  direction).  URRDF  * URIDF  cut  0.  For  a Lambertian  surface.  BRRDF  e 
Const,  occurs.) 

i(|v  Component  of  she  URIDF  contributed  by  reflection  from  the  rough  air.matetul  interface. 

Preceding  page  blank 
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Symbol  Definition 

w Lambertian  component  of  the  BR1DF. 

Pj_  Lambertian  component  of  the  directional-hemispherical  reflectance,  which  is  the  ratio  of 

the  total  radiant  power  [W]  reflected  into  all  directions  to  the  radiant  power  incident  [W| 
(from  a chosen  direction). 

<?in  The  minimum  value  found  during  optimization  for  a.  the  rms  of  the  difference  between  the 

theoretical  and  experimental  values  of  the  BR1DF  divided  by  the  theoretical  value. 

INTRODUCTION 

In  this  report,  we  wish  to  present  some  contributions  to  the  understanding  of  the  reflection  of 
light  from  rough  surfaces.*  For  simplicity,  we  consider  only  rough  surfaces  that  are  flat  on  a macro- 
scopic levei  but  rough  on  a microscopic  level.  The  microscopic  irregularities  are  assumed  to  have  ran- 
dom shape  and  distribution  over  the  surface  and  no  preferred  orientation  m any  direction  along  the 
plane  of  the  macrosurface  (a  uniform,  symmetric,  rough  surface).  For  perspective,  the  following  brief 
summary  of  the  current  understanding  of  light  reflection  from  such  surfaces  is  presented 

Lambert's  Lew 

The  firsl  theory  of  light  scattering  by  a rough  surface.  Lambert's  Law,  can  be  derived  by  con- 
sidering the  incident  light  ray  to  he  “completely  randomized”  by  the  rough  surface,  tty  completely 
randomized  is  meant  a purely  hypothetical  situation  m which  an  incident  ray  emerges  Iroitt  the  surface 
(unallected  by  the  air-material  interface)  after  last  being  scattered  isotropically  and  with  equal  likelihood 
from  every  point  tit  the  semiflnite  space  below  the  surface  plane  U tuber  f's  Law  states  that  the  reflected 

radmtee  (reference  l>  (W  st‘l  (proj.  cm*  )'1  j (or  the  famtltar  property  of  visual  brightness)  is  invariant 

with  observation  direction,  and,  for  constant  incident  ir radiance,  |Wcm  *]  is  invariant  with  incidence 
direction.**  Also.  Lambert's  Law  states  that  the  reflected  light  o unpolart/ed  and  independent  of  in- 
cident polarization  Howcvet.  many  surfaces  do  not  even  vaguely  follow  Lambert's  law  ’die  best 
diffusing  surfaces  known  deviate  a few-  percent  (tefetences  2 and  ,1)  even  at  zeto  zenith  angles  of  in- 
cidence and  reflection,  and  differ  grossly  at  zenith  angles  peater  than  -15  degrees  Tbete  3tc  two 
mechanisms  that  tend  in  randomize  the  incident  light  and  thus  tend  to  produce  Lambertian  cflevtton 
multiple  reflection  from  the  irregularities  of  the  rough,  au  inatertal  interface  and  multiple  scattering 
from  inhomogencmes  w-fthm  the  material  below  the  interface  The  multiple  reflection  mechanism  t.\ 
considered  secondary’,  since  ii  o found  (reference  4)  fo  contribute  only  a few  percent  to  tltc  reflection 
from  iih4'  dielectric  surfaces.  Most  of  this  reflection  may  he  from  mils  double  reflection  and  thus  not 
randomized  compel  civ. 


“Tint  uerfc  was  supported  fey  the  Naval  Air  Ssitcmt  f'cuitmand  unde:  AlKT  \Sh>  \ mini  ,’)S  *ti|  in.1t  t-ont  \u 
Weapons  System)  InvsSlipatioos.  and  A J*$>.V'225/WdO-?44*{!t . Sutfa.c  Tatyct  (.rudarcc  idiotic  (.uulatwr 
ApplK-aitom) 

**As  an  aid  In  the  trader,  rhe  ur.its  of  cash  letro  or  tprantils  still  folio*  n m hiackcit  Area  «mt»  preceded  hs  "pmj  ” 
refer  to  the  area  of  a ssrtfaec  projected  onto  a plane  perpendicular  to  the  p*np*eaU««  dues  non  »f  she  lieht  Olhci 
1 area  unUs  tc(e»  to  an  area  upon  a uirfime 
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The  multiple  scattering  is  the  primary  mechanism  for  randomizing  the  incident  light  and  acts  as 
follows.  Light  entering  the  surface  material  is  somewhat  randomly  scattered  by  refraction  at  the  rough, 
air/material  interface.  Then  the  light  is  randomly  scattered  one  or  more  times  by  inhomogeneities  within 
the  material.  Finally,  it  is  again  somewhat  randomly  scattered  by  refraction  upon  exiting  the  interface. 
(A  similar  mechanism  has  been  evaluated  numerically  in  reference  5.)  Thus,  three  or  more  random 
scatterings  must  occur,  and  three  should  greatly  randomize  the  light.  However,  light  is  completely  ran- 
domized only  after  an  infinite  number  of  random  scatterings;  thus,  some  non-Lambertian  light  could 
emerge.  Also,  it  has  been  shown  (reference  6)  that  the  light  becomes  partially  polarized  (in  contradic- 
tion to  Lambert's  Law)  upon  exiting  the  rough,  air/material  interface  and  that  such  polarization  can  be 
quite  pronounced  at  large  zenith  angles  of  reflection.  Thus  the  light  can  be  partially  de-randomized 
upon  exiting  the  interface. 

Lommd-Saeliper  Law 

The  above  described  multiple-scattering  can  be  divided  into  scattering  at  the  interface  and  scattering 
below  the  surface.  The  below-surface  scattering  should  contribute  the  most  to  the  randomization  of  the 
multiple-scattered  light.  Below-surface  scattering  has  been  treated  extensively  (see  references  7.  8,  and  9 
for  summaries).  The  Lommel-Seeliger  Law  was  derived  by  considering  only  single,  isotropic  scattering 
by  below-surface  scatterers  and  below-surface  attenuation  before  and  after  scattering.  Kven  though  light 
icflected  in  this  way  is  not  at  all  randomized,  it  was  found  to  give  only  a slight  deviation  from  Lambert’s 
Law.  The  Lommel-Seeliger  Law  was  extended  by  including  n 'nisotropic  scattering  (various  “phase  func- 
tions" for  the  scatterers)  and  double  scattering  Chandrasekhar  (reference  10)  succeeded  in  including  ah 
orders  of  multiple  scattering  below  the  surface.  His  results  for  some  phase  functions,  especially  the 
isot topic,  did  not  deviate  much  from  Lambert's  Law  except  at  latge  zenith  jngtes  of  incidence  and  re- 
flection. 

Shadowing  of  parts  of  a Lambertian  surface  by  the  surface  irregularities  was  also  considered  (ref- 
erences 7 and  1 1).  Such  shadowing  lends  to  produce  a peak  reflection  m the  direction  of  die  source, 
but  is  not  usually  a large  effect. 

Thus  the  above  mechanisms  may  not  produce  perfectly  Lambertian  lellectnm.  hut  the  on  per  fee 
turns  cited  cannot  fully  account  for  the  large  deviations  from  Lambert's  l aw  found  for  many  surfaces 

Hqjh-  and  low  RefUxtwc*  Surfaces 

Hark,  or  low-reflectance,  dielectric  vurtaces  are  found  (o  deviate  the  most  Itiom  1 amheri's  bawl 
For  dielectrics,  the  reflection  from  the  air  matenal  interlace  is  low  ta  few  petcentl  and  sancs  little  tot 
various  materials  The  rest  of  the  incident  light  (~**5  percent  I passes  into  the  surface  material,  where  it. 
and  thus  the  reflectance,  o strongly  affected  bv  the  absorption  coefficient,  which  vanes  greatly  with  the 
material.  If  the  coefficient  is  high . the  below-surface  reflectance  can  be  zero,  leaving  only  the  interface 
reflectance,  which  is  always  low  Thus,  for  Sow-reflcs lance  surfaces,  reflection  from  the  interface  is 
much  nv>re  significant-  Also,  flits  reflection  is  found  (reference  4j  to  be  almost  entirely  single  reflection, 
which  is  not  at  all  randomized  and  thus  might  deviate  grossly  from  I amber  t'v  law 

Bougucz  Theory 

A theory  for  single  t-fectmn  by  i rough.  iiMmieital  interface  was  futi  presented  by  Bouguci 
(reference  ID  about  Pw.i  The  interface  is  modeled  as  csuisisnng  of  randomly  oiicnted.  optically  flat 
facets,  each  of  which  behaves  Hkc  a small  plane  nmtoi.  reflecting  a pisrtmu  of  the  lighi  me  idem  on  it 
as  determined  bv  Fresnel's  equations  foi  specular  reflection  at  a flat  dielectric  interface  For  cidlimated 
incident  light,  those  facets  of  ivnty  one  orrematn  n can  reflex t light  into  a given  direction  Thus,  all  tight 
in  that  directum  has  been  reflected  in  exactly  the  same  was . so  n »» not  at  ail  randomized . and  the  re- 
flection by  the  facet  mechanism  is  not  at  all  l. amber  nan 
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The  Bouguer  facet  theory  has  been  found  to  explain  much  of  the  non-Lambertian  reflectance  and 
has  undergone  considerable  development.  (1)  The  derivation  of  the  basic  theory  has  been  refined  (re- 
ferences 13  through  19).  (2)  Several  models  (references  13,  17,  18,  and  20  through  22)  have  been  trieu 
for  the  “microarea  distribution  function."  This  is  the  function  that  gives  'he  relative  number  of  facets 
oriented  in  any  given  direction,  or,  more  precisely,  the  relative  total  facet  surface  area  per  unit  solid 
angle  of  surface  normals  pointed  in  any  given  direction.  This  function  governs  the  directional  distribu- 
tion of  the  scattered  light.  If  the  function  is  uniform  with  respect  to  the  direction  of  the  facet  normals, 
a diffuse-like  distribution  occurs.  If  most  of  the  facets  lie  nearly  parallel  to  the  plane  of  the  surface,  a 
nearly  specular  disttibution  occurs.  (J)  Another  development  of  the  facet  theory  is  the  inclusion  of 
shadowing  (references  II,  14,  17,  19,  and  23).  that  is,  some  facets  cannot  contribute  to  the  reflection 
because  they  are  in  the  shadow  i/  neighboring  irregularities  The  shadowing  theory  results  m a "geo- 
metrical attenuation  factor,"  the  reflection  in  any  direction  being  reduced  by  a complicated  function  ul 
only  the  direction  of  incidence  and  reflection.  Shadowing  has  been  found  (reference  19)  to  be  the  only 
mechanism  able  to  explain  the  sharp  cutoff  of  light  near  grazing  reflection.  (4)  Another  and  recent  de- 
velopment is  the  extension  of  facet  theory  (reference  1 9)  to  cover  reflection  outside  the  plane  f inci- 
dence. However,  there  has  been  little  comparison  with  measurements  outside  of  the  plane  of  incidence. 

Other  Mechanisms 

There  are  some  other  mechanisms  possibly  contributing  to  rough -suf lace  reflection.  There  may  he 
randomly  oriented  facet  -like  areas,  “rough  facets."  each  of  which  reflects  as  a little  flat  l atnhetinin  re- 
flector (reference  16).  resulting  on  more  light  being  soltored  at  higher  zenith  angles  ol  letleetion  than 
from  the  ideal  large  flat  l.amhertian  suiljce.  Hits  might  be  the  erne  if  the  below -mu  face  scattering 
originated  mostly  from  a depth  th it  was  much  smaller  than  the  typical  size  of  the  facets 

There  may  be  reflection  from  a pl.u,,,.  1 wet  boundary  to  the  rough-sutface  material,  .is,  for  ex 
ample,  for  painted  flat  surfaces.  This  is  found  ireierence  24)  to  explain  some  reflection  phenomena 

There  ntav  be  Bougucthke  facet  reflection  Itom  below -surface  cracks,  flits  t;  pc  of  reflection 
produces  effects  differing  from  those  ol  surface  Ucei  reflection  and  helps  to  explain  some  phenomena 

(rel'creme  25) 

Diffraction  Mechanism* 


The  flouguet  mechanism  and  the  aho.se  mechanisms  are  based  on  t,r\  optics  and  thus  feipme  that 
the  size  of  the  surface  irregularities.  and  all  then  ra-Jit  of  curvature  be  much  greater  than  the  wavelength 
ol  the  light  Itelcreneev  2ts  and  2')  Tim.  ot  course.  may  tun  alwass  be  tine,  so  some  ot  the  tefleefton 
may  be  explainable  soils  by  s»ax<  optics  Beckmann  ami  Sptz/iehim>  Deference  2b)  give  the  basic  theory 
ol  the  application  of  diffraction  theory  to  reflcvtn'n  lonn  ihe  tough  interlace  of  a randomly  tough  sur- 
face. Vans-  physical)'  important  tevtihv  ate  oh  tamed  tn  ch«ed  fosm.  hot  (hex  are  complex  and  limited, 
for  example  infinite  sene'  often  result,  lew  polarization  phenomena  can  he  treated,  and  sha  towing  has 
in  be  ignored  One  scry  stgttillc.ini  result  of  (he  effects  of  wave  optics  is  the  Raxlotgh  Criterion  tretcr. 
cnees  14  and  2**)  This  criterion  implies  i|u(  at  chough  angles  to  grazing  incidence  on  a rough 
surface.  which  at  a.tt  other  angles  ol  incidence  teflecis  according  to  ray  optics,  there  xx til  appear  a 
specti!*;  'eiHvtroii  va**scd  bx  diffraction  The  large'  the  megulantiex  and  then  radii  of  curvattiie.  the 
vniallei  the  rebutted  grazing  angle  xvill  hv-.  but  the  diffraction-caused  specular  reflection  will  iiifMnatch 
appear  Relevant  to  tax  reflectr.sn  at  the  ait 'material  interlace.  Beckmann  obtains  a closes!  fi'tni  solution 
m the  large  itrcgularuy  limit  for  a surface  characterized  by  a <«auxxian  shsirrbuirot;  >.<  surface  heights  and 
an  auto  coriclatts'n  length  along  Ihe  surface  Sudi  a haracicnzalnm  merely  imply s aii«tlv„*t  micmarca 
distributism  function 
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On  most  veal  surfaces,  the  microstructure  probably  includes  some  parts  whose  reflection  can  be 
approximated  by  ray  optics  and  other  parts  whose  reflection  must  be  treated  by  physical  optics.  One 
would  suspect  that  usually  the  broad  facet-like  areas,  being  large  and  flat,  could  be  treated  by  ray  optics, 
but  that  the  corners  and  edges  of  these  areas,  being  small  and  highly  curved,  must  be  treated  by  physical 
optics.  Much  recent  work  (references  30  through  37)  has  been  done  on  diffraction  from  edges  and 
corners,  resulting  in  some  simple  expressions  for  the  reflected  light  distribution  and  some  ways  for 
experimentally  separating  th^  diffracted  light  from  the  ray-reflected  light.  The  diffraction-scattered 
light  was  predicted  to  be  qui'e  specular. 


Rigorous  Solutions  of  Maxwell's  Equations 


When  the  wavelength  becomes  large  enough  relative  to  the  size  of  the  scattering  irregularity,  dif- 
fraction Ibimulae  such  as  those  of  Beckmann  and  Spizzichino  (reference  28).  based  as  they  are  on  a 
scaler  wave  equation,  predict  results  that  deviate  significantly  from  the  results  obtained  from  a rigorous 
application  of  Maxwell's  equations.  This  topic  has  been  the  subject  of  much  effort  (and  some  contro- 
versy) in  the  IELL  Journal  of  Antennae  and  Propagation  during  1971,  1972.  and  1973. 


A Generalization 


Although  deviations  in  some  directions  and  for  many  surfaces  are  large,  it  appears  (references  14, 
17.  20.  38,  and  39)  that  most  of  the  reflection  in  most  directions  for  most  rough  sut faces  is  closely 
approximated  by  a combination  of  Lambertian  reflection  from  below  the  interface  and  simple  (no  shad- 
owing) Bouguer  reflection  from  interface  facets.  This  is  reasonable  theoretically  because:  diffraction 
phenomena  are  generally  limited  to  a narrow  region  near  the  specular  direction;  shadowing  produces 
extremely  large  effects  only  near  tow  elevation  angles  of  reflection:  multiple-surface  scattering  is  found 
to  be  small : below-stirface  particle  scattering  tends  to  be  nearly  Lambertian;  reflection  from  a hclow- 
surface  plane  should  usually  be  small  because  of  attenuation  in  the  surface  material;  mirror  reflections 
from  helow -surface  facets  will  be  somewhat  randomized  from  entering  and  exiting  the  rough,  air/ material 
interface:  and  the  "rough-facet"  reflection  cannot  give  really  radical  departures  from  Lambertian 
reflection. 


New  Model  for  Ray-Tiieory  Reflection 


1 Ins  report  presents  a new  model  for  ray -theory  reflection  from  the  rough,  air/material  inteiface. 
Instead  of  randomly  oriented  flat  microarcas  (facets) , the  surface  structure  is  here  modeled  as  an  en- 
semble of  randomly  oriented,  randomly  curved  microareas.  Since  for  practical  purposes  any  curved 
area  can  be  broken  down  into  infinitesimal  facets,  the  two  models  are  identical  (physically)  and  thus 
must  give  the  same  results.  However,  this  new  deriv.uion  leads  to  some  new  and  useful  results;  that  is, 
a new  visualization  of  the  surface  structure  (an  “average"  irregularity)  is  conceived  and  proved  valid;  a 
new.  mine  aecuiate.  and  useful  function  representing  the  surface  structure  is  discovered  (resulting  from 
an  ellipsoid  of  revolution  as  the  "average"  irregularity);  and  a way  to  normalize  the  resulting  reflection 
equations  to  give  the  absolute,  instead  of  just  the  relative,  reflectance  distribution  is  found. 


This  model  is  then  comMue1’  with  some  Lambertian  reflection  and  tested  on  data  from  a partic- 
ular application  of  imeiesi.  The  application  of  interest  is  to  laser  target  designation  systems.  A laser 
illuminates  a point  on  a taiget  and  a icceiver,  on  a piece  of  ordnance  or  on  a fire  control  system,  de- 
tects the  direction  of  the  illuminated  point  and  directs  the  ordnance.  The  geometry  for  this  application 
consists  of  incidence  from  any  angle  of  a polarized  .ollimuied  beam  onto  a flat  rough  surface;  observa- 
tion fiom  a distance  such  that  the  illuminated  ana  is.  for  practical  purposes,  a point;  and  observation  of 
the  reflected  light  from  any  direction,  hot!,  within  and  without  the  plane  of  incidence  (we  believe  that 
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this  is  the  first  time  a theory  of  ray  reflection  front  the  interface  is  extensively  tested  outside  the  plane 
of  incidence).  Tite  model  is  limited  to  s>mply  a combination  of  ray  reflection  front  the  rough,  air/ 
material  interface  and  Lambertian  reflection  be  ause  this  combination  has  been  found  (references  14, 
17,  20,  38,  and  391  to  give  most  of  the  reflection  in  most  directions  front  most  surfaces,  and  because, 
for  the  application  to  target  designation  systems,  a model  is  needed  that  gives  only  reasonable  accuracy 
for  a wide  variety  of  common  rough  surfaces  and  uses  s few  surface-related  parameters  as  possible. 

The  comparison  to  the  data  was  good.  Discrepancies  can  generally  be  explained  by  gross  deviations  of 
some  of  the  measured  surfaces  from  the  assumed  surface,  such  as  the  existence  of  two  surface  materials 
instead  of  just  one.  ar  dm  existence  of  a significant  directional  dependence  of  the  surface  structure. 
The  validity  of  the  n->  ait,  .-.ion  is  verified,  and  the  usefulness  of  the  concept  of  an  ellipsoid  of  revo- 
lution as  the  average  surface  irregularity  is  demonstrated. 


Because  of  the  geometry  described  above,  the  usual  and  most  generally  useful  quantity  for  express- 
ing the  reflectance  distribution  is  poorly  applicable,  so  a less  generally  used  but  more  applicable  quantity  is 

adopted  Nicodemus (reference  40)  calls  the  unusual  quantity  the  BRRDF,  f , , [W  sr  "* (proj.  cm* )-1  / W cm"  | 
(bidirectional  refiected-radiance-distribution  function),  or  for  short,  the  BRDF.  f (bidirectional  reflectance- 


distribution  function).  It  is  defined  as  the  radiance  [W  sr_l(proj.  cm")  * ] reflected  into  a given  direction 


per  unit  irradiancc  (W  cm""]  incident  on  the  surface  from  another  given  direction.  This  quantity  involves 
radiance,  which  is  very  useful  because  it  has  certain  invariance  properties  (reference  41)  upon  transfor- 
mation through  optical  systems  and  upon  reflection  from  surface  to  surface.  However,  irradiance  is  a 
density  of  light  incident  on  a surface,  and  for  our  geometry,  a reflectance-distribution  function  containing 
such  a density  is  not  the  most  useful  quantity.  Much  better  for  our  purposes  is  the  BRIDF.  f j,  (bidirec- 
tional reflccted-intensity-distribution  function)  (reference  42).  It  is  defined  as  the  radiant  intensity  [W  sr-1 1 
reflected  into  any  specified  direction  from  an  illuminated  "point"  per  unit  radiant  flux  (W | in  the  incident 
beam.  To  find  the  more  generally  useful  BRRDF  (or  BRDF).  simply  divide  the  BRIDF  by  the  cosine  of 
the  zenith  (or  sine  of  the  elevation)  angle  of  relied  ion. 


Organization  of  Report 


This  report  is  organized  as  follows.  For  the  derivation  of  the  interface  reflectance  model,  first  the 
reflection  from  a single  curved  microarea  is  derived,  then  an  ensemble  of  microareas  is  formed,  resulting 
in  a reflectance-distribution  function  in  terms  of  a surface  structure  function,  and  the  concept  of  an 
"average"  irregularity  is  presented  and  validated.  The  model  is  put  m terms  of  the  laboratory  coordinates, 
and  the  Fresnel  reflectance  of  the  surface  material  is  found.  Nest,  a normalization  procedure  is  imple- 
mented to  make  tite  re  ll  ee  l a nee  d i s I nln  1 1 ion  function  an  absolute  quantity,  and  an  ellipsoid  of  revolution 
is  used  as  a one  parameter  model  for  the  average  irregularity.  For  the  comparison  to  measurements,  some 
Lambertian  reflection  is  added,  and  the  resulting  four-parameter  tellecdon  model  is  applied  to  reflectance- 
distribution  function  measurements  both  within  and  without  the  nominal  plane  of  incidence  on  a variety 
of  rough  surfaces. 


INTERFACE  REFLECTION  MODEL 


Reflection  From  a Single  Microarea 


Consider,  as  shown  in  figure  I . a very  small  curved  area  AAv  on  alt  irregularity  ot  «h-  tough-surface 

mierosiruciure.  This  area  is  chosen  small  enough  that  we  can  approximate  the  curvatures  in  all  di.ections 
on  it  bv  sections  of  circles.  (The  curvaitues  of  the  circle  sections  may  he  different  in  different  directions 


across  the  area.)  Also,  the  area  is  chosen  small  enough  that  the  angles  spanned  by  these  circle  sections  are 
very  small.  This  smallness  allows  the  specification  of  the  orientation  of  be  a single  suitace  normal  it. 


The  direction  ol  n is  given  (rela  ivc  to  the  planar  macrosurlace)  hv  a zenith  angle  « ard  an  azimuth  angle 
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Figiioj  1.  Reflection  From  a Curved  Microarea  of  the  Rough-Surfaco  Air/Motorial  Interface. 
A small.  optically  smooth,  slightly  curved  area  AA,  with  a normal  n in*ercepts  the  light  (rum 
an  area  AA  of  the  Incident  beam  and  diverges  it  Into  Aiu  story Jiam  in  a direction  r such  that 

ft  A A A 

r A it  « a v f. 


i defined  as  shown  in  figure  t.  Let  fhe  incident  light  be  a narrow  collimated  beam  originating  front  a 

A 

source  in  the  direction  /.  and  let  the  incident  light  hmv  a radiant  flux  (reference  I)  of  *l*|  (Wj  spread  (not 
necessarily  uniformly)  over  the  beam's  cross-sectional  area  / At  the  point  in  A,  where  AA,  is  located, 
let  the  incident  normal  irradiance*  (flux  pet  unit  area  projected  (H»rpond tcular  to  the  beam)  be  Ln  |W  (proj. 
cm-1)-'  j . 

Let  AA,  be  optically  smooth  and  let  its  minimum  dimension  and  its  minimum  radius  of  curvature  be 
much  greater  than  the  wavelength  of  the  incident  light,  Then  (references  2t»  and  27)  ray  theory  applies,  so 


•Imdtanee  |W  cm"1 1 h the  tins  incident  upon  a linn  area  of  a surface.  In  an  unvhoojtiniJ  incident  beaut,  the  irradiance 
varies  with  the  orientation  of  the  unlace  tanple  of  incidence!  and  is  thus  not  a constant  of  the  beam.  "Normal  irradiance" 

| W (prop  cm-’r1 1 is  the  irradiance  upon  a surface  oriented  normal  to  the  incident  beam  and  thus  is  a constant  of  beam. 


;-'f'v":  *»»wc, . -*.  ■ ■ : ’£. . i^wAs’"- 

■■*—  ■■■■•■"• ■"  -•■•'  ’’  ' <<x  /*  ’ -v  * _ • '■}'  y'd&  - ■ 


AAS  reflec's  as  a section  of  a curved  mirror  surface:  that  is,  of  the  incident  light  intercepted  by  AAS,  a por- 
tion &{s,n,k)  (determined  by  the  Fresnel  equations:  n and  k are  the  refraction  and  absorption  indices. 


respectively)  is  reflected  into  a small  cone  (represented  for  simplicity  by  the  rectangular  cross  section  cone 


in  figure  ! ) of  solid  angle  Aw  (sr]  in  a direction  r determined  by  the  mirror  reflection  equation  r x it  = ft  x /. 


Th"  refle  ted  l ^it  is  not,  strictly  speaking,  in  the  form  of  a cone.  It  would  be  in  the  form  of 
a true  cone  ii  the  curvature  of  AAs  did  not  vary  with  the  direction  on  AAS,  since  AAS  would  then  be 

a spherical  mirror,  and  the  reflected  light  would  radiate  from  a definite  point  image  (real  or  virtual  de- 
pending on  the  sign  of  the  curvature).  The  vertex  of  this  cone  would  lie  along  the  reflection  direction 
at  a distance.  ,>efore  or  after  AA^,  determined  by  the  magnitude  and  sign  of  the  curvature  of  AA_.. 

But  since  the  curvature  in  general  does  vary  wi'h  direction  on  AAs,  the  image  point  is  smeared.  But 
the  ra>s  refected  from  AAs  (or  their  backward  extensions  in  the  case  of  a virtual  image)  all  pa«s 

through  a football-shaped  volume  of  space  approximately  bounded  in  one  direction  by  two  points  and 
in  the  other  directions  by  the  circumference  of  an  area  whose  si/e  is  directly  related  to  the  size  of 
AAs.  The  two  bounding  points  are  the  image  points  for  spherical  mirrors  of  curvatures  equal  to  the 
maximum  and  minimum  curvatures  found  on  AA^.  Thus,  except  for  the  realistically  impossible  cir- 
cumstance of  a perfectly  zero  maximum  ,-r  minimum  curvature,  the  dimensions  of  this  volume  arc  finite. 
Now,  let  the  observe;  be  at  infinity.  To  him.  thi-  volume,  me  smeared  image,  is  effectively  a point,  and 
thus  the  reflected  light  is  effectively  a true  cone.  Also,  this  condition  makes  the  position  of  the  image 
point  irrelevant  (for  convenience,  let  the  cone  vertex  lie  on  AAs). 


We  can  now  quantify  the  tight  -ciioctcd  from  AAv  as  a radiant  intensity  (ladiant  Dux  per  stetadian 
radiating  from  a point  source),  since  the  light  relieved  from  AA^  originates  effectively  fium  a point. 
This  reflected  radiant  intensity  |\V  sr  1 ] is  derived  io..n  the  curvatures  o.'  AAs  as  follows.  Of  the  in- 
cident beam  cross  sectional  atea.  the  area  intercepted  by  AAv  is  AAs  cos  s |cnr| . and  thus  AAv  telle;  :r 
a radiant  flux  of  <>($)  fen  AAS  cos  s (\V|  into  the  small  solid  angle  Aw  |si| . giving  the  reflected  rad'am 
intensity  I s 6?(s)  |*n  AAs  cos  s/Aw  |W  st*1 1 , ...cuing  AAx  be  small  enough  that  it  subtends  only  a small 
angle  in  any  direction  allows  the  curved  lines  A/(l  and  A J(  (see  figure  2i  to  be  very  nearly  stiaight  lines. 
Thus,  the  approximation  AA,.  = A lu  A S holds,  where  AJ,  is  the  line  segment  formed  on  A £ by  the 

intersection  of  AAv  with  a plane  parallel  to  the  macrosurface  and  A I is  the  tine  s.gment  formed  on 
, A 

by  the  intersection  ol  wuh  the  plane  containing  the  mucrosurhue  tunmal  k ami  the  JuWur- 

face  normal  n.  Since  these  line  segments  ate  sections  of  circle.,  they  can  he  given  h>  A/,  = u A/I.  and 

A Ja  * 9 |An|.  where  oi  and  o(|  (defined  positive)  ate  tne  tadii  of  curvature  of  AAx.  and  A/  and  Au  are 

the  arcs  subtended  by  A/;  rrJ  £J  . respective!).  Strictly,  u is  not  a radios  of  cutvulure  of  AA^  because 

it  is  not  perpendicular  to  AA^,  hot  this  does  not  affect  the  derivation.  The  reflected  solid  angle  can  be 

given  by  Aw  * cos  0 1A 0 Ad/ 1 (sr) . Substituting  the  above  three  expressions  into  I gives 


1 - 


<{(s.n.k)  l:n  cos  s , AaA, 


COS  {’ 


A0  Ad/ 


ou  ot  (W  sr’1 1 


III 


Since  o and  / are  functions  of  0.  d/.  and  jl-  the  quantity  lAu  A/’Ad  AvI  can  be  given  by  the  Jacobian 
determinant  J(ay.:d.d»l  (in  the  limit  as  the  incremental  angles  apptouch  zeto). 
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Figuitl  2.  Specification  of  (he  Shape  ol  a Curved  Micthjte*  in  Term*  ol  Two  OMhojonal  Curved  Line* 
A ihjhtly  curved  elemeni  ol  *urt»ce  area  AA^  may  be  miertected  by  a plane  parallel  to  (he  macro- 

turloce  »o  lorm  dw  hue  with  a curvature  radio*  ol  a?.  and  AA^  may  he  nKerteciad  by  (he  plane 
parallel  (o  n and  £ to  form  (he  line  A/  with  a curvature  rodiui  ol  a 

<1  ft 


Since  «.  /.  s,  and  J ate  functions  ot  only  the  lahotalotv  coordinates  (ti,  (>.  and  ami  n and  k arc 
parameter*  ot‘  only  the  surlace  moienal.  the  only  quantities  in  equation  ( I)  that  depend  on  the  size  or 
shape  of  AAx  arc  «a  and  <>r  Thus,  Ok  surface  geometry  appear*  in  equation  i I ) only  as  the  curvatures 

of  AAx.  Also,  the  sine  of  AAx  drws  nor  enter  and  may  be  disregarded. 

Note  that  the  derivation  is  still  valid  if  either  one  or  both  ol  the  lines  St  and  A / have  a curvature 

/ 13 

opposite  that  shown  in  llgute  This  fact  lets  equation  (I)  apply  to  a curved  area  of  any  shape  convex, 
concave.  <>t  sjddlc.  Also,  since  radii  of  curvature  are  defined  as  positive  and  the  Jacobian  determinant  is 
used  only  in  an  absolute  value.  I remains  unchanged  sf  the  curvatures  of  either  A/  or  A / ot  both  are 
reversed. 

£nieenbic  of  Microamt  sod  the  Surface -Structure  Function 

Consider  lire  following  model  for  the  surface  mtcrostructure  l.et  the  microsurface  be  continuously 
curved  and  randomly  undu'atmg.  similar  to  lull  or  uumiium  topography . Let  it  l>c  optically  smooth  and 


let  all  curvature  radii  and  irregularity  sizes  be  much  larger  than  the  incident  light  wavelength.  Let  the 
AAs  be  any  small  area  on  this  undulating  surface.  Under  these  conditions,  equation  (1)  gives  the  reflec- 
tion from  each  small  area.  Equation  (I ) may  be  applied  to  a discontinuous  microstructure  if  one  disre- 
gards diffraction  effects  at  corners  and  edges.  This  allows  application  of  equation  (1)  to  a randomly 
scratched  surface,  such  as  ground  glass  or  toughened  metal,  and  to  a “globule  pile"-like  structure,  such  as 
some  diffuse  spray-painted  surfaces.  Also,  equation  (1)  may  be  applied  to  a surface  containing  flat  facets 
if  one  makes  their  curvatuie  radii  not  quite  infinite.  A very  wide  variety  of  microstructures  is  thus  in- 
cluded. Finally,  let  the  statictieal  properties  be  uniform  with  position  on  the  surface  and  independent  of 
direction  along  the  surface  (no  wood-like  grain  or  scratches  with  a preferred  direction). 

For  given  directions  of  illumination  and  observation,  an  observer,  if  his  visual  resolution  were  sharp 
enough  to  resolve  the  surface  microstructure,  would  see  the  reflected  light  originating  only  from  various 
points  scattered  around  the  microstructure.  These  are  points  at  which  the  microsurface  normal  happens 
to  be  in  the  proper  direction  for  reflection  to  the  observer.  Since  the  illuminated  macrosurface  area  is 
small  and  the  receiver  is  at  infinity,  the  relative  locations  of  the  reflecting  points  are  irrelevant,  and  all 
the  reflected  light  may  be  treated  as  if  it  originated  from  a single  point.  Thus  we  may  speak  of  a re- 
flected total  radiant  intensity  lT  {W  sr-1 1 equal  to  the  sum  of  the  radiant  intensities  contributed  by  the 
individual  reflecting  points.  Normals  at  the  reflecting  points  must  all  be  parallel,  so  the  reflection  from 

each  point  is  given  by  equation  <1 ) with  the  same  values  of  C.  C».  a.  / and  s.  If  for  a particular  combina- 
tion of  incidence  and  reflection  directions,  a given  incident  beam  illuminates  N reflecting  points,  the  total 
radiant  intensity  reflected  in  the  specified  direction  is  given  by  the  ensemble 


(if(s.n.k)  h cos  s 
cos  0 


(W  sr-’  | . 


( 1 a) 


This  equation  nia>  he  simplified  by  replacing  the  normal  n radiance  t |W(proj.  cm*  )_l  | incident  on  each 
reflecting  point  by  an  avetage  normal  u radiance  given  by  En  ( = «l»j  ‘Aj  |W  cm”'| , where  «l>(  |W|  is  the  radiant 
flux  of  the  incident  beam  and  A(  is  the  ci oss-sectional  aiea  of  the  incident  beam  This  does  not  require  iliat 

flic  incident  beam  be  unifoun.  but  only  ( 1 1 that  a statistically  large  number  ol  reflecting  points  he  illumina 
ted  and  (1)  that  the  incident  beam  he  uniform  enough  so  that  no  large  portion  of  the  incident  flux  lulls 
on  only  a lew  reflecting  points.  The  total  radiant  intensity  tints  becomes 


’l  a 


4\(s.n.kl  <|>.  cos  s 
cos  0 


■fe)l  it  -- 

' i"  I 


<:> 


It  all  dependence  on  the  surface  material  and  laboratory  geometry  can  be  separated  into  one  factor 
and  all  dependence  on  the  surface  structure  into  a second  factor  ibis  second  factor  would  be  a surface- 
structure  function.  Notice  that  all  dependence  rtf  the  surface  structure  is  localized  in  the  summation 
factor  of  equation  (2);  however  this  factor  contains  a dependence  on  the  laboratory  geometry  , namely. 

N varies  with  incidence  angle  ;)  This  can  be  removed  as  follows:  For  normal  incidence  (incidence  angle 
0 a 0)  and  a chosen  direction  of 'reflection, consider  the  set  of  N reflecting  points.  Their  curvatuie  radii 

products  (ofli  o . t 3 1.2 N’nl  compromise  the  summation  factor  in  equation  (2).  As  the  incidence 

angle  increases,  this  same  set  of  points  reflects  into  a new  direction.  The  new  radiant  intensity  is  given 
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by  equation  (2)  with  different  values  of  6,  y,  a,  J.  etc.,  but  with  the  same  set  of  oQi  azj  in  the  summa- 
tion factor.  However,  as  the  angle  of  incidence,  p,  increases,  the  illuminated  macrosurface  area  increases, 
as  1/cos  p,  and  some  points  in  the  additional  area  also  contribute  to  the  reflection.  So  the  summation 
factor  for  the  new  radiant  intensity  must  contain  some  additional  aaiazj  quantities.  Assuming  that  the 

surface  is  statistically  uniform,  these  extra  quantities  will  be  statistically  the  same  as  the  original  set,  and 
the  summation  factor  is  merely  increased  by  1/cos  p.  Thus  equation  (2)  becomes 


tf?(s,n Jc)  4>j  cos  s 

I = J 

cos  p cos  8 


Iwsr~’i 


and  the  incidence  angle  dependence  has  been  separated  from  the  summation  factor.  However,  the  summa- 
tion factor  still  retains  a dependence  on  tne  laboratory  geometry,  namely,  Nn  is  proportional  to  the  cross- 

sectional  area  Aj  of  the  incident  beam.  This  is  easily  compensated  for  by  including  the  1/Aj  factor  with 
the  summation  factor. 


Thus,  the  quantity 


D(u)  = 


II 


°ai  azi 


[dimensionless] 


is  a surface-structure  function,  since  it  contains  all  the  dependence  on  the  surface  structure  and  no  de- 
pendence on  anything  else.  In  general,  this  quantity  is  a function  of  both  the  coordinates,  a and  z.  of 
he  microarea  nc  mal,  but  we  are  considering  only  surfaces  whose  statistical  properties  are  directionally 
iniform,  and  for  such  surfaces  D would  have  no  z-dependence.  In  the  Bouguer  facet  theory  and  its 
refinements,  the  surface-structure  dependence  is  incorporated  as  a microarea-distribution  function,  which 
is  the  relative  amount  of  microarea  oriented  in  a given  direction  or  the  probability  density  of  a facet  nor- 
mal to  be  in  a given  direction.  Like  D,  this  function  contains  all  the  dependence  on  the  surface  structure 
and  no  dependence  on  anything  else.  Since  the  two  reflectance  theories  are  identical  physically,  this 
proves  that  D and  the  microarea  distribution  function  must  be  the  same,  or  at  least  proportional. 

Equation  (3)  can  be  put  in  terms  of  a reflectance-distribution  function  by  dividing  by  the  incident 
flux  <!>,  [W| . This  gives  the  PRfDF  (reference  42)  ff!s  [sr-1 1 . Incorporating  equation  (4)  in 

addition  gives 


fT[s(p£,\p)  = <R(sjv.k) 


cos  p cos  0 


m 


[W  sr"  /W| . 


(Division  by  sin  8 gives  the  BRRDE  (reference  42)  [W  sr"1  (proj.  cm2)_,/W  cm"2] ). 
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Average  Surface  Irregularity 

It  might  prove  useful  if  there  could  exist  a single  optically-sinooth  curved  surface  that  would  reflect 
light  in  the  same  distribution  as  that  reflected  by  the  ensemble  of  all  the  curved  microareas  comprising 
the  rough  surface.  This  would  mean  that  a randomly  irregular  surface  could  be  treated  as  if  it  consisted 
of  a large  number  of  small  identical  “average"  irregularities  or  that  a randomly  irregular  surface  could  be 
treated  as  if  it  were  a single  large  curved  surface  (for  a uniformly  intense  incident  beam  of  the  same  total 
flux).  Such  an  average  irregularity  would  have  to  be  a surface  of  revolution  about  the  macrosurface 
normal  because  of  the  assumed  directional  independence  of  the  rough  surface  statistics. 

Mathematically,  the  question  of  the  possible  existence  of  such  an  average  irregularity  may  be  stated 
as  follows.  Let  pa(a)  and  pz(a)  be  the  radii  of  curvature  for  the  average  irregularity  (defined  as  were, 

respectively,  oa  and  oz.  for  a curved  microarea),  and  let  C be  a constant.  There  must  exist  a surface  of 
revolution  whose  D(a),  given  by 

D(a)  = C,oa(a)pz(a)  (dimensionless] , (6) 

is  the  same  as  the  D(a)  resulting  from  equation  (4).  It  may  be.  but  is  not  obvious,  that  for  any  physically 
realistic  functional  form  of  10(a),  given  by  equation  (4),  there  exists  a surface  of  revolution  which,  by 
equation  (6),  can  give  this  same  functional  form  of  D(a).  It  will  be  proved  that  this  is  indeed  true  and, 
furthermore,  that  the  shape  of  the  surface  of  revolution  can  be  very  greatly  restricted  and  still  be  general 
enough  to  be  able  to  give  any  physically  realistic  functional  form  of  D(a). 

The  following  restrictions  may  be  applied  to  D(a)  to  make  it  conform  to  physical  reality. 

1.  D(a)  must  exist  for  all  values  of  a at  least  between  0 and  tr/2.  Any  realistic  rough  surface  will 
have  some  microareas  at  any  given  value  of  a between  0 and  n/2. 

2.  D(a)  must  be  finite  for  all  values  of  a,  since  no  real  surface  would  contain  any  perfectly  flat  or 
perfectly  cylindrical  microarea.  (In  wave  optics,  all  fiat,  rough  surfaces  produce  a perfectly  specular 
reflection,  of  some  magnitude.  This  produces  a delta  function  in  the  reflectance-distribution  function 
and  thus  an  infinite  D(a)  at  one  point.  However,  we  are  treating  only  ray  optical  reflection.) 

3.  D(a)  is  positive  for  all  values  of  a.  This  follows  from  equation  (4),  considering  that  a.,  and  aui 
were  defined  such  that  they  were  always  positive. 

4.  D(a)  must  obviously  be  sii.Je-valued.  No  physical  quantity  can  have  more  than  one  value  for  itself 
at  the  same  point, 

5.  D(a)  is  continuous.  Physical  quantities  never  have  perfect  discontinuities. 

First,  it  will  be  proved  that  there  exists  a surface  of  revolution  that  can  give  not  only  the  functional 
forms  of  D(a)  that  might  result  from  equation  (4).  but  also  any  physically  realistic  functional  form  of 
D(a).  Let  h(x)  be  a curve  that  when  revolved  about  x = 0 gives  the  surface  of  revolution  (see  figure  3). 
The  p7  curvature  radius  is  given  by  p?  = lx|.  The  p„  curvature  radius  is  the  radius  of  curvature  of  h(x), 

so  it  is  given  by  pa  = (1  + (h')2 ! 3/2/|h"|  (panics  indicate  derivatives  with  respect  to  x).  The  angle  a is 
equal  to  the  negative  slope  angle  of  h(x);  that  is,  a - tan_,(-h').  Substituting  these  expressions  into 
equation  (6)  gives 
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Figure  3.  The  Concept  of  an  "Average"  Surface  Irregularity.  The  "average"  surface 
irregularity  is  generaterf  by  rotating  the  curve  h{ x)  about  x = 0.  It  is  an  optically  smooth 
surface  of  revolution  of  a shape  such  that  it  gives  the  same  distribution  of  reflected 
light  (when  irradiated  by  a uniform,  well-collimated  beam)  as  that  given  by  the  actual 
rough-surface  microstructuro  (when  irradiated  by  the  actual,  nonuniform,  well-collimated 
beam).  The  shape  of  h(x)  may  be  restricted  to  the  general  shape  shown  and  still  be 
general  enough  to  represent  any  microstructure.  This  general  shape  consists  of  the 
slope  being  zero  ut  x = 0 and  infinity  at  h - 0 and  the  curve  between  having  no  inflection 
points  or  straight-line  sections. 

Substitution  of  p(x)  = h'(x)  gives 

lp'l/!x|  = C(1  + p2)-V2/D|tan'1(-p>] . 
which  is  equivalent  to 

p’  = Cx  (1  + pV/2/D[tan_l(-p)|.p7x  >0. 
p’  = -Cx  (I  + p2 >3/2/D |tan”*(-p)I • p'/x  < 0. 


Equations  (7b)  and  (7c)  are  of  the  form  p'  = f(x,p).  “A  differential  equation  of  this  form  has  a solution 
p = p(x)  through  every  point  (x  * x0,  p = p())  with  a neighborhood  throughout  which  lfx.p)  is  continu- 
ous." Since  the  tfx.p)  of  each  of  the  above  equations  is  everywhere  continuous,  all  points  (x0.  p())  have 
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continuous  neighborhoods,  and  a solution  p = p(x)  to  each  equation  exists  for  all  values  of  x.  Since  both 
equations  (7b)  and  (7c)  have  solutions  for  all  values  of  x,  equation  (7a)  must  have  a solution  p = p(x) 
for  all  values  of  x.  Since 

h(x)  = /p(x)dx  + C, 

(Cj  is  an  arbitrary  constant),  and  since  p(x)  exists  for  all  values  of  x,  then  h(x)  must  exist.  Thus,  there 

must  exist  a solution  h = h(x)  to  equation  (7)  for  any  physically  realistic  D(a);  that  is,  there  must  exist 
a surface  of  revolution  capable  of  representing  any  physically  realistic  rough  surface. 

With  this  established,  it  will  now  be  shown  that  not  only  is  there  an  h(x)  for  every  D(a);  but  that 
also  some  limits  to  the  form  of  h(x)  can  be  established  such  that  h(x)  is  still  capable  of  giving  any  form 
of  D(a);  that  is,  because  of  the  nature  of  equation  (7),  certain  options  are  allowed,  the  choice  of  one  of 
which  will  not  limit  the  ability  of  h(x)  to  give  any  D(a).  Two  options  are  the  arbitrary  choices  of  two 
boundary  conditions  for  a second-order  differential  equation:  one  condition  on  h'(x)  and  one  on  h(x). 
This  is  allowed  because  a solution  to  a differential  equation  is  obviously  s1  ill  a solution  when  restricted 
by  a boundary  condition,  and  the  set  of  solutions  h(x)  for  all  forms  of  D(a)  are  obviously  still  solutions 
for  their  respective  forms  of  D(a)  even  though  the  same  boundary  condition  is  applied  to  every  h(x).  A 
third  option  is  the  choice  of  a plus  or  minus  sign  for  h”(x).  This  is  allowed  because  only  its  absolute 
value  appears  in  equation  (7).  Thus  we  have  established  some  ways  of  limiting  the  shape  of  the  surface 
of  revolution. 

Thus  far  we  have  found  two  types  of  ways  of  limiting  the  functional  form  of  h(\)  while  not  limiting 
its  ability  to  give  any  physically  realistic  fin  nr  of  D(a).  These  were  three  arbitrary  options  of  restrictions 
on  h(x)  resulting  from  the  nature  of  equation  (7)  and  five  restrictions  on  D(a)  resulting  from  physical 
reality.  These  will  be  applied  in  the  following  way  to  determine  a very  limited  type  of  curve  for  h(x). 
but  a type  of  curve  that  is  still  versatile  enough  to  give  any  physically  realistic  form  of  D(a).  First,  let 
the  boundary  condition  on  h'(x)  be  h'(0)  = 0.  This  makes  the  slope  of  h(x)  zero  at  x = 0.  Second,  use 
the  option  on  the  sign  of  h"(x)  to  make  h”(0)  negative.  This  forces  the  slope  of  h(x)  to  begin  to  de- 
crease as  x begins  to  increase  (from  zero).  Third,  the  h(x)  cannot  have  inflection  points.  It  is  did.  the 
same  value  of  a would  occur  at  different  parts  of  the  curve.  Since  D could  have  different  values  at  these 
two  points,  there  would  be  two  values  of  D for  the  one  value  of  a.  Thus  D would  be  multiple-valued, 
and  restriction  (4)  would  be  violated.  Fourth,  the  h(x)  cannot  be  a straight  line  at  any  point.  If  it  were, 
Pa  would  be  infinite  at  that  point  and  thus,  by  equation  (6),  D would  be  infinite,  and  restriction  (2) 

would  be  violated.  Thus  far,  the  slope  has  been  forced  to  be  level  (zero)  at  x = 0 and  start  to  slope  down 
(decrease)  as  x starts  to  increase.  Since  h(x)  can  have  no  inflection  points  or  straight-line  segments,  the 
slope  must  continue  to  decrease  with  increasing  x.  Thus  'he  slope  must  become  vertical  (-°°)  at  a finite 
value  of  x (say.  x„).  [f  the  slope  only  approached  vertical  as  x approached  infinity,  h(x)  would  approach 
a straight  line  and  D would  approach  infinity,  and  restriction  (2)  wo„,d  be  violated  (sec  the  appendix  for 
a rigorous  verification  of  this  contention).  Last,  let  the  boundary  condition  on  h(x)  be  h(x())  » 0.  This 
forces  the  vertically  sloping  part  of  h(x)  to  lie  on  the  x-axis. 

In  summary,  a limited  functional  form  of  h(x)  capable  of  giving  any  physically  realistic  functional 
form  of  D(u)  resembles  the  curve  in  figure  A.  This  curve  is  level  at  x = 0.  decreases  with  increasing  x 
without  inflection  points  or  straight  sections,  and  becomes  vertical  as  it  crosses  the  x-axis.  Or.  in  terms 
of  a.  the  curve  is  such  that  a ~ 0 at  x = 0,  and  a increases  without  stopping  or  turning  back  witli  increas- 
ing x to  become  a = at  h - 0. 

It  is  instructive  to  be  able  to  visualize  the  random  topography  of  a rough  surface  as  being  equivalent 
to  one  large  surface  of  revolution  of  very  restricted  shape,  but  what  substantive  contribution  diws  this 
make?  It  does  not  put  any  limits  on  the  functional  form  of  the  surface-structure  function  D(a),  as  we 
have  just  proved.  But  it  may  give  insight  into  the  discovery  of  better  models  to  represent  the  surface 
structure,  indeed,  as  is  shown  later,  the  most  obvious  choice  of  a surface  of  revolution,  an  ellipsoid  of 
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revolution,  gives  a surface-structure  function  that  is  better  than  the  existing  ones.  Perhaps  other  shapes 
for  the  surface  of  revolution  will  give  even  better  results.  Also,  maybe  ones  involving  two  or  more  param- 
eters will  give  very  good  higher  order  approximations.  To  derive  the  D(a)  for  any  surface  of  revolution, 
the  procedure  used  later  for  the  ellipsoid  can  be  followed. 

Model  in  Terms  of  Measurable  Variables 

For  the  result,  equation  (5),  to  be  useful,  the  unmeasurable  variables  a.  z,  and  s and  the  Jacobian  J 
must  be  found  in  terms  of  the  measurable  variables  ft,  6,  and  ift.  The  angles  a and  z can  be  found  from 
the  reflection  equation 


A a A a 

r x n = n x / 


(8) 


since  this  is  the  condition  that  forces  n(ay.)  to  be  in  a determined  direction  for  light  from  a source  in 
the  direction  Hft)  to  be  reflected  in  the  direction  r(fl.iji).  These  three  vectors  are  expressed  as  follows: 

A A A 

/ = j sin  ft  + k cos  ft,  (*)) 

A A A A 

n = i sin  t sin  a + j cos  / sin  a + k cos  a.  ( 10) 

f = f sin  t li  cos  0 + j cos  0 cos  0 + Is  sin  0 . (II) 

Substituting  these  equations  into  equation  (H)  gives  the  following: 

cos  / sin  a cos  ft  - cos  u sm  ft  - cos  >>  cos  0 cos  u - sin  0 cos  i.  sin  a,  (12) 

-sin  / sin  a cos  ft  = sin  0 sin  /.  sin  a - sin  y>  cos  0 cos  a.  ^ (13) 

sm  /.  sin  d = sin  y cos  0 cos  / - cos  0 cos  0 sin  /.  (14) 

The  a(d.O.O)  and  /(d.d.y)  can  be  found  by  solving  equations  ( 13)  and  ( 14),  respectively.  Then 

tan  a 1 Mil  y cos  fl'[sin  / (sin  0 * cos  d)|.  (15) 


tan  / = sin  sos  0; (sin  ft  * cos  y cos  0).  (16) 

Hquution  (l<»)  gives  zdi.O.C'l.  and  if  this  is  used  m equation  (13).  equation  (IS)  gives  tUd.^.y).  The 

a ,• 

angle  s can  be  found  by  substituting  equations  (‘J>  and  (10)  into  the  relation  /•  n ■ cos  s: 

cos  s 13  cos  / sin  u sin  (i  ♦ cor  a cos  d (l?l 

Since  o(d.d.O)  and  /(d.O.yl  are  given  hy  equations  (13)  and  <IM,  tins  gives  std.O.y).* 


♦Similar  expressions  are  derived  l>\  torrenxe  and  Spaimsv  irelctoiue  tM|  tn  iismparc  ihcec  and  nilwi  results.  the 
lotluwine  vomet'inn  heiv.  een  our  and  I licit  symbolisms  wa»  found 

t S R li  « v ' «. 

I4S  „ till  - tl  0 4 » v.  ’ -*t 

When  equations  US).  < I A I.  and  1 17 1 were  converted  in  ihe  ViS  symbolism  and  enmpared  in  the  cottespondinc  es- 
preishms  nl'  T&S.  they  wete  I'nund  in  he  very  dmimdar  in  form.  Ilnwevet.  numerically  lliev  .impaled  idem  walls 
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With  the  use  of  equations  (15)  and  (i6).  the  Jacobian  determinant  J can  be  evaluated: 
-cos2z  cos2a  sin  0 cos20  (1  + sin  0 cos  0 + cos  0 cos  t//  sin  (3) 


J 


fe)- 


sin  z (sin  0 + cos  0)2  (sin  0 + cos  t J/  cos  0)2 


This,  however,  was  found  numerically  to  be  equal  to 


J 


(5) 


-cos  0 


4 cos  s sin  a 

Although  it  was  difficult,  this  equivalence  was  also  verified  analytically. 


(18) 


(19) 


«$k~ 

*v 


Micro vm  Reflection  Coefficient 

The  fi(s,njc)  in  equation  (S)  can  be  found  from  the  Fresnel  equations  (with  a complex  index  of 
refraction).  The  reflectances  <5^  and  di|j  for  incident  light  linearly  polarized  perpendicular  and  parallel  to 

the  local  plane  of  incidence,  respectively,  arc  given  (references  43  and.  in  somewhat  different  form, 
reference  44)  by  the  following: 


sin*(s  - s,)  ♦ sinlry 

4i|(sjrJk)  » . CO) 

«n*(s  ♦ s,>  ♦ smh’x 


fi|j(sjr Jt)  e 


cos*(s  ♦ Sj)  ♦ smh'x 
ccvs2(s  - M ♦ smhJx 


Cl) 


where 


cos  Sj  a b/sinhx . 


with 


b * k sin  s/(n2  ♦ k2). 
c * 1 - un2*/in2  ♦ k2). 
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and  n and  k.  res  -ctively,  the  real  and  imaginary  components  of  the  complex  index  of  refraction  of  the 
surface  material.  The  k is  sometimes  called  the  absorption  index  of  the  surface  material.  Of  academic 
interest  are  s1  and  x.  which  come  from  Snell’s  Law  for  a complex  index  of  refraction. 


sin  s/sinfs,  + ix)  = n - ik 


At  s = 0.  equations  (20)  and  (21 ) are  indeterminate.  However,  the  solutions  are  found  (reference  45)  to 
be 

n2  + k3  + l-2n 

4(|(0.njc)  = «|j(0jiJs)  = — . <--) 

n~  + k~  + I + 2n 


The  above  expressions  give  6?  for  light  polarized  relative  to  the  local  plane  of  incidence;  however, 
for  laboratory  measurements  6\  must  be  determined  lor  light  polarited  relative  to  the  macroscopic  plane 
of  incidence,  since  this  is  the  usual  way  that  the  incident  polarization  is  specified.  The  angle  (set ween 
the  macroscopic  and  local  planes  of  incidence  gi  is  the  angle  between  their  respective  normals,  t and  / x it 
(see  figure  4).  so  <p  is  given  by 

A {'  A, 

cos  0 - i 'lx  it/ sin  s. 

Substituting  equations  (*>)  and  (10)  into  the  above  equation  gives 

cost?  * (sin  d cos  u - cos  (i  cos  / snt  uVstn  s.  (25 1 

for  which  everything  on  the  tight-hand  side  has  been  found  The  orientation  of  the  incident  polarization 
plane  is  specified  by  the  angle  y.  the  angle  that  the  incident  polarization  plane  makes  with  the  macroscopic 
plane  of  incidence  The  angle  y - £ is  then  the  angle  that  the  modem  polarization  plane  makes  with  the 
local  plane  of  incidence,  hot  calculation  of  the  reflectance  of  the  incident  Itglii.  ns  polarization  mum  be 

resolved  Into  components  perpendicular  and  parallel  to  the  local  plane  of  nteidence.  A part  sttt*(y  - P) 
of  the  incident  light  ts  polarized  perpendicular  to  the  local  plane  of  incidence  lot  which  the  reflectance 

is  <1^1  s).  and  a pari  ct*s‘(y  - d)  of  the  incident  light  ts  polarized  parallel  to  the  local  plane  nf  incidence 
for  which  the  reflectance  if  *s yt > > Thus  the  reflectance  of  light  potan/ed  at  an  attgfc  y from  rhe  macro- 
wopic  plane  of  inc*  4*nce  is 

^s.y)  = *lj(s)  unMy  - $)  * em’ty  - C'  <24f 

Prom  this,  the  refiectatvce  of  light  polatt/ed  parallel  and  petpendtcular  t»  the  macroscopic  plane  of  in- 
cidence is  given  by 

<Jv(s)  ■ *l(s.y  8 0)  = 4fjUI  sni'i?  ♦ 


4<Mfs»  * «(s.y  * *)2)  = lljts)  ctn*0  ♦ 4{(jts)  stn*0. 


respcctivdy 
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Ktgkjr*  4 R<r*tt<*n  Between  *h»  Local  »«<j  M»c»<nu<l»ce  Plane  o»  ItKntamtn  The  vector*  / amt  e he 

J * 

in  'he  m*t»Qui«fae«!  plan*  at  motSOKe.  »n<l  I an<t  h he  >n  the  loco  til  ane  o'  mtulente  The  an#ie 
between  I hew  plane*  n ecjiiai  n>  o lihe  an^le  txuwwtn  the.*  no«n*i*l 

live  feflpcUtwe  of  unpolatwed  of  ctreulatly  pohtn/ed  light . o simply  the  average  oi  4i,ut 
anti  since  hi  th  unpolatoed  ami  atcuUtlv  poiatved  light  tcvolve  mio  emu  I components  petpett 

dicular  and  parallel  to  the  Uveal  plane  ot  incidence 


l^tsl  ♦ «iHU)| 


The  above  equations  tcveal  no  inhumation  on  the  polati/atton  lotm  ot  the  reflcctcsl  light  If  the  te 
fleeted  polatwation  form  tv  dented,  (torn  Holl  (tcfetcncc  44).  the  phase  shift  -Jiffeicme  between  h'. 

and  41[|  u given,  ami  horn  Jvhutehff  (tefctence  4<>)  the  Muelct  mains,  ( son  taming  41, . 4?^.  f>.  and  t>)  hat 
the  polatitattvm  transformation  on  relleetum  can  be  derived 
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Normalization 

Substituting  equations  (19)  and  (6)  into  equation  (5)  gives 


frls(0£,£>  = <R(s.74i,k) 


C p_(a)  p fa) 


4 cos  (i 


[W  sr'Vw], 


sin  a 


(25) 


As  follows,  the  constant  C can  be  found  so  as  to  make  the  reflected  intensity-distribution  function  fr)s 

an  absolute  quantity.  Since  'be  average  irregularity  is  a surface  of  revolution,  the  region  around  a = 0 
can  be  treated  as  a spherical  mirror  surface  of  radius  pQ(0).  For  normal  incidence  and  reflection,  it  is 

easy  to  derive  the  BRIDF  for  a spherical  mirror.  This,  of  course,  must  be  equal  to  the  BR1DF  of  our 
model  for  aouital  incidence  and  reflection  frJs  (0°5)Q°,  0°),  and  C can  be  found  from  this  equality. 

The  focal  length  of  a spherical  mirror  is  equal  to  half  the  radius,  or  —p^O),  and  collimated 
incident  light  reflects  into  a cone  with  the  focal  point  as  its  vertex.  A flux  of  Chj/Aj)  AA  [W|  is  in- 
cident on  a small  area  AA  on  the  top  of  a hemisphere.  If  the  single  large  surface  of  revolution  is  to 
replace  all  the  illuminated  areas  of  the  real  surface,  it  must  have  uniformly  intense  light  incident  on  all 

parts  of  it  and  have  no  light  missing  it.  Therefore,  the  area  of  its  base  npj(n/2)  [cm2]  must  equal  the 

incident  beam  cross-sectional  area  A,,  and  thus  <t>,  AA/(np^(n/2))  [W]  is  the  flux  incident  on  AA.  Upon 

1 t 

reflection,  this  flux  is  attenuated  by  <R(0)  and  is  diverged  into  AA/[—  Pp(0)]“  stcradians.  Then  (R(0) 
multiplied  by  the  above  incident  flux  and  divided  by  .his  number  of  stcradians  gives  the  reflected  radiant 
intensity  [W  sr_l  ].  Additionally  dividing  this  by  <h[  [WJ  gives  the  BRIDF  for  the  spherical  mirror: 


frls  (0°90°.  0°) 


<K(0)p;(0)/(4rrp^(rr/:)l  [W  Sr-'/W| 


(26) 


For  toe  rough-surface  model  at  normal  incidence  and  reflection,  (3  = 0,0  = n/2.  a - 0,  s = 0,  and 
P7(0)  = 0 occur,  and  equation  (25),  the  BRIDF  becomes 


f . (0h90°,  0°)  = -7C  <R(0)  p (0)  lim  (p  (a)/sin  a)  [Wsr'/Wj. 

r<s  ^ a a-0  z 


(27) 


Equating  equations  (26)  and  (27)  and  solv.ng  for  the  normalization  constant  C gives 


(.’  = p„(0)/jrrp;  (r,/2)  lim  (p  (a)'sin  a)| . 

1 a—O  '■ 


(28) 
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Substituting  this  into  equation  (25)  gives  the  normalized  or  absolute  form  for  the  BRIDF.* 


1 Pj°)  Ptt(«)  P,(«) 

- fl(s,y,n,k)  — [W  sr'  /WJ  (29) 

7f  cos  p p2  lint  (p  (c)/sin  a)  sin  a 
1 a-0  i 

i 

i 

Ellipsoid  of  Revolution  Average  Surface  Irregularity  Approximation  | 

\ 

I 

It  is  not  practicable  to  make  direct  measurements  of  the  oaj  and  azi  values  or  of  pQ(a)  and  pz(a). 

The  best  that  one  can  do  is  to  choose  a reasonable  model  irregularity.  A sphere  is  one  possibility,  but 
it  contains  no  parameter  that  one  can  vary  to  change  the  characteristics  of  the  surface  structure.  An 

ellipsoid  of  revolution,  h = e(l  -x2)>,!,  has  one  such  parameter  e,  which  is  the  ratio  of  the  length  of  the 
axis-rotated-about  to  the  length  of  the  axis-rotated.  As  e decreases,  the  model  irregularity  becomes  ; 

flatter  and  reflects  more  light  in  the  specular  direction.  This  is  a very  useful  property,  since  the  wide  |_ 

variety  of  real  surfaces  contains  a continuous  distribution  of  diffuse-to-highly-specular  surfaces.  The 
f .(0 #,'/')  for  the  ellipsoid  of  revolution  is  found  by  evaluating  equation  (29)  as  follows.  The  derivative  ' 

of  h(x)  for  the  ellipsoid  is  ■ 

h'  = -ex(l  - x2)"^  = - tana. 

Solving  for  x gives  I 

x = tar.  a (tan2a  + e2)~'/>  = pt(a) 

I 

The  second  derivative  of  h(x)  is 

h"  = -e(l-x2)"3^2  = -e[l  - tan2a/(tan2a  + e2)]-3^2. 

j 

The  pa(a)  is  the  radius  of  curvatuie  of  h(x)  and  is  thus  given  by  pa(a)  = [1  + (h')2]3^2/lh''l,  Substitut- 
ing h'  and  h''  this  equation  gives 

pja)  - e2  (( 1 + tan2a)/(e2  + tan2a)]3^2. 


The  other  components  of  equation  (29)  are  found  to  be: 


p (0)  = 1/e ; p,(n/2)  = 1 ; and  lint  [p,(a)/sin  a!  = l/e. 
s 1 a-*o  z 

Substituting  the  above  expressions  into  equation  (29)  gives  the  BRIDF  (reference  42)  for  the  ellipsoid  of 
revolu  tion  -ave  rage  i rregu  I art  ty . 


♦One  could  extend  this  model  to  include  shadowing  by  simply  multiplying  equation  (29)  by  Torrence  and  Sparrow's 
shadowing  factor  (symbolism  converted  by  use  of  the  table  in  the  footnote  following  equation  (17)).  This  would 
not  affect  the  normalization,  since  their  shadowing  factor  is  normalized  to  equal  ono  when  no  shadowing  occurs. 
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WW)  = 


fl(s,7,n,k) 


4ir  cos  0 (e2  co$2a  + sjn2a^2 


(W  sr“*/W] 


(Division  by  sin  6 gives  the  BRRDF  [W  sr"1  (proj.  cm2)-,/W  cm-2) .) 

Thus,  with  equation  (30)  we  have  a simple  model  giving  the  absolute  BRIDF  contributed  by  ray 
reflection  from  the  irregular  microsurface  * Only  three  parameters  are  involved:  the  refractive  index  n 
of  the  surface  material,  the  absorption  index  k of  the  surface  material,  and  the  axis  ratio  e of  the  ellipsoid 
of  revolution-average  surface  irregularity  representing  the  surface  structure.  AH' other  quantities  in  equa- 
tion (30)  have  been  found  in  terms  of  the  laboratory  coordinates  0.9.  and  \ jj. 


COMPARISON  TO  MEASUREMENTS 

Earlier,  we  speculated  that  the  concept  of  an  average  surface  of  revolution  representing  the  surface 
structure  might  give  insight  into  the  choice  of  a better  surface-structure  function.  The  first  choice,  the 
ellipsoid  of  revolution,  resulted  in  the  surface-structure  function  given  by  the  last  term  in  equation  (30). 
Normalized  to  D(a  = 9)  = 1 , it  is 


D(a)  = e4/(e' 


cos2a  + sin2a)2 . 


(31) 


Since  this  surface-structure  function  was  shown  to  be  equivalent  to  the  microarea  distribution  function 
used  as  the  surface-structure  function  in  the  Bouguer  facet  theory,  this  function  can  be  compared  directly 
to  the  others.  Rense  (reference  21)  found  a way  to  determine  the  microarea  distribution  function  from 
the  reflectance-distribution  function  measurements.  Figure  5 presents  his  data  for  one  rough  surface, 
along  with  plots  of  best  fits  for  our  structure  function  and  three  microarea  distribution  functions  found 
in  the  literature.  The  much  closer  fit  to  the  data  shows  our  function  to  be  a significant  improvement. 

To  test  (2-e  accuracy  of  the  model,  to  illustrate  the  value  of  the  concepts  derived,  and  to  test  the 
applicability  of  'he  model  to  laser  target  designation  systems,  the  model  was  combined  as  follows  with 
some  Lambertiar  reflation  and  compared  to  BRIDF  measurements  made  on  a variety  of  rough  surfaces. 

The  BRIDF  contributed  by  the  Lambertian  reflection  is  simply  given  by 


frlL(3^^)  = sin  0 pl(0-2n)lTr  [sr~* ) . 


(32) 


where  pL(0j2n)  (Lambertian  directional-hemispherical  reflectance  (references  42  and  47)(  is  the  portion 

of  the  incident  power  scattered  by  Lambertian  reflection.  This  quantity  might  vary  significantly  with 
incidence  angle,  but  for  lack  of  any  relationship,  it  will  be  assumed  constant.  This  is  added  to  equation 
(30)  to  give  the  combined  BRIDF: 


t fr,(/JjM>  - fr,s(W)+  (sr-'|. 


(33) 


♦Equation  (30)  was  found  to  compare  to  Torrence  and  Sparrow's  result.  To  compare  the  two  results,  our  symbolism 
was  transformed,  the  two  surface  structure  functions  were  made  equal  to  one,  ’I  AS’s  shadowing  factor  was  removed, 
and  the  BRIDF  was  transformed  to  T&S’s  ratio  of  BRRDFs,  This  completed,  the  results  were  the  same. 
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RELATIVE  MICROAREA  DISTRIBUTION  FUNCTION 


O Data  taken  from  Rente  (reference  21) 

A 7 O O O 

' Plot  of  our  function  e He  cos  a + sin  a)  (e  ■=  0.40) 

2 2 2 2 

*••*•••••••  Plot  of  function  e /(e  cos  a + sin  a)  (a  = 0.25),  originated  by  Berry  (reference  13) 

_4  2 2 2 

Plot  of  function  (cos  a)exp(-A  tan  a)  (A  = 7.62)  derived  from  Beckmann's  (reference 
n i • 28)  result  for  a surface  characterized  by  a Gaussian  distribution  of  surface  heights  and  an 

autocorrelation  length  in  the  ray  optics  approximation 


-2  2 2 2 

Plot  of  function  (cos  a)exp(-A  tan  a)  (A  = 6.93)  used  by  Sirohi  (reference  22). 
originated  by  Berry  (reference  13! 

2 2 2 2 

Plot  of  an  approximation  exp(-A  a ) (A  ■ 0.0021  /degrees’  ) of used  by  Rense 

(reference  21)  for  small  values  of  a 


a,  ZENITH  OF  MICROAREA  NORMAL  (DEGREES) 
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Figure  5.  Comparison  of  Best  Fits  of  Some  Surface-Structure  Functions  (Equivalent  to  the  Microarea 

or  Facet  Distribution  Function). 
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A trial  and  error  optimization  procedure  was  applied  to  relative  BR1DF  measurements  made  on  a 
number  of  surfaces.  This  optimization  procedure  consisted  of  a method  of  varying  the  four  parameters 
[e,  n,  k,  and  pL(d^ir)|  of  the  model  (equation  33)  until  the  minimum  value  am  was  found  for  a,  which 

is  the  root  mean  square  normalized  deviation  of  the  experimental  values  from  the  theoretical  values. 


,M 

l i=! 


where 


it  '-)/(l  4 


Xj  is  measured  relative  value  of  relative  BRIDF;  fr]j  is  the  theoretical  value;  and  M is  the  number  of 

measurements  on  a given  surface.  The  set  of  measurements  for  each  surface  consisted  of  the  following. 
For  each  combination  of  values  of  6 of  15,  30.  45.  60.  75.  and  lH)  degrees  and  of  (J  of  0,  15.  30.  45. 

60,  and  70  degrees,  several  values  of  4>  were  chosen,  and  the  relative  BRIDF  was  measured  for  collimated 
6,328-Angstrom  continuous  wave  laser  light  polarized  perpendicular  and  parallel  to  the  plane  of  incidence. 
This  resulted  in  about  400  measurements  for  each  surface.  This  number  is  statistically  large,  so  the  loca- 
tion of  om  in  the  space  of  the  four  parameters  should  not  be  sensitive  to  small  random  experimental 
errors. 

Table  1 gives  a and  the  corresponding  (optimum)  values  of  the  parameters  for  eh  surface. 

Figure  6 compares  a small  sample  of  the  measurements  to  the  model.  A brief  description  of  each  surface 
and  a discussion  of  the  results  is  given  in  the  following  paragraphs. 

The  3M  Black  Velvet  paint  is  a very  dull  antireflection  paint  for  optical  instruments.  The  optimiza- 
tion resulted  in  two  very  nearly  equal  minima  of  a.  The  associated  two  sets  of  values  for  the  parameters 
of  the  model  give  nearly  identical  reflection  patterns;  therefore,  reflection  patterns  are  plotted,  in  figure 
5,  for  only  the  first  set  of  parameters  listed  in  table  1.  The  model  applied  quite  well.  The  curves  com- 
pared very  well  with  the  measured  points  in  shape  and  magnitude.  Only  the  dependence  on  incident 

Tibi*  1.  Optimum  Valuta  of  tha  Modal  Paramatart 


Paranwttr 

3M  Black  Valvat  Paint 

Oirtiad  Oliva 
Drab  Paint 

Cament 

Plywood 

Orass 

anr> 

0.290* 

0,293* 

0,345 

0.200 

0.186 

0328 

0.020 

0.022 

0.06 

0.33 

0.40 

0.5 

e 

0.89 

0.88 

0.24 

0.68 

0,51 

1.6 

n 

1.60 

1.62 

1.0 

1.00 

0.9 

0.7 

k 

i 

0.00 

0.3 

0,4 

0.9 

1.9 

1 3 

•Two  minima  (or  o were  found. 
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*«fft  “ 0403 
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0.03  0.03  0.01 


to  tfl  - 0.63  et  * - 180°  (x  art  ffJ  » 1.88) 


to  ff(  ■ j80^  ^ 


0.07  003  006  004  0.03  002  0.01  0 0.01  0.02 


0.30  0.25  930  0.16  0.10  0.06  0 0.06  0.10 


x at  ifl  - 08; 


0.7  08  08  0.4  0.3  0.2  0.1  0 0.1  0.2 


0.6  08  0.4  0.3  0.2  0.1  0 0.1  03 


by  the  radius  vector  from  the  origin  to  the  point  or  curve,  and  C'  is  given  by  the  polar  angle.  (The  BRRDF 
(reference  421,  reflected  radiance  per  incident  irradiance,  is  given  bv  fr( /sin  ().]  Values  for  0 and  0 for  each 

column  of  graphs  are  given  at  the  top.  The  symbols  x and  © indicate  measured  values  of  f j,  and  the 
continuous  solid  and  dashed  curves  give  the  theoretical  plots  of  f | for  the  incident  light  polarized 
respectively  perpendicular  and  parallel  to  the  plane  of  incidence. 
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polarization  gave  a consistent  error.  However,  this  error  was  small  for  most  points.  Also,  from  the  ab- 
solute values  of  BR1DF  given  by  the  model,  the  directional-hemispherical  reflectance  (sum  of  flux  scattered 
in  all  directions  per  incident  flux  from  one  direction)  at  normal  incidence  can  be  estimated  and  was  found 
to  be  of  the  order  of  magnitude  of  4 percent.  This  is  consistent  with  the  measured  directional-hemispher- 
ical reflectance  value  of  2.5  percent.  Also,  under  a microscope,  the  surface  appeared  as  an  irregular  pile 
of  little  spheres,  which  is  consistent  with  the  optimum  ellipsoid  axis  ratio  of  e = 0.89,  which  is  very 
nearly  that  of  a sphere  (e  = 1).  The  difference  could  be  caused  by  the  sphere  sagging  slightly  into  ellipses 
when  they  were  still  liquid. 

The  dirtied  olive  drab  paint  was  a sample  of  a slightly  glossy  dark  olive  drab  paint.  It  was  soiled  by 
placing  dry  soil  on  its  slightly  dampened  surface  and  shaking  off  all  that  would  not  stick.  The  model 
applied  rather  poorly.  The  shapes  of  the  curves  and  the  dependence  on  incident  polarization  were  only 
roughly  similar,  and  the  magnitude  was  systematically  erroneous  for  many  curves.  However,  the  directional- 
hemispherical  reflectance  at  normal  incidence,  estimated  from  the  absolute  values  of  BRIDF  given  by  the 
model,  is  of  the  order  of  5 percent,  which  is  consistent  with  the  measured  value  of  4 percent.  Also,  the 
ellipsoid  axis  ratio  e is  very  small,  as  it  should  be  for  a glossy  surface:  that  is,  a glossy  surface  should 
have  a large  portion  of  its  microsurface  area  with  normals  near  the  microsurface  normal,  as  does  a very 
flat  ellipsoid  of  revolution.  Much  of  the  poor  performance  of  the  model  might  be  because  there  are  two 
types  of  reflecting  surfaces:  the  paint  and  the  soil.  The  optimization  assumed  only  one. 

The  cement  surface  consisted  of  common  structural  concrete  cement  troweled  flat  and  allowed  to 
harden  in  a horizontal  position.  The  surface  appeared  to  be  slightly  varied  in  shades  of  light  grey.  The 
model  applied  very  well.  There  were  no  systematic  errors.  The  small  random  errors  could  easily  be 
attributed  to  the  varied  nature  of  the  surface,  since  the  illuminated  area  of  the  surface  was  not  necessarily 
the  same  for  each  measurement.  This  is  consistent  with  the  fact  that  no  such  random  errors  occurred 
during  the  parts  of  the  measurement  program  in  which  the  illuminated  area  remained  the  same.  The 
illuminated  area  remained  the  same  for  only  cases  in  which  either  the  incident  polarization  was  varied  or 
when  \p  was  varied  at  0 =»  0.  The  directional-hemispherical  reflectance  at  normal  incidence,  estimated 
from  the  absolute  values  of  BRIDF  given  by  the  model,  is  of  the  order  of  50  percent,  which  is  consistent 
with  measured  values  ranging  from  30  to  50  percent  for  various  samples  of  concrete. 

The  plywood  surface  was  clean  and  unaged.  It  had  an  irregular  grain  whose  width  was  about  equal 
to  the  diameter  of  the  laser  beam.  A large  portion  of  the  grain  was  oriented  nearly  parallel  to  the  fibers 
of  the  wood.  Throughout  the  measurements,  the  plane  of  incidence  remained  parallel  to  the  fibers.  The 
modei  applied  fairly  well.  The  directional-hemispherical  reflectance  at  normal  incidence,  estimated  from 
the  absolute  values  of  BRIDF  given  by  the  model,  is  of  the  order  of  70  percent,  which  is  reasonable, 
since  the  surface  appears  to  be  highly  reflecting.  Also,  lire  lowest  o of  airy  surface  occurred,  and  the 

grainy  nature  of  the  surface  could  easily  account  for  the  little  remaining  deviation  if  the  deviation  were 
random.  However,  much  of  the  deviation  causing  c was  systematic.  The  dependence  on  incident 

polarization  did  not  compare  very  well,  and  the  experimental  points  deviated  in  the  form  of  a bulge 
whose  maximum  appeared  to  shift  from  v = 90°  to  0 - 180°  as  (3  varied  from  0 to  90  degrees.  This 
type  of  systematic  error  would  result  from  this  particular  directional  dependence  of  this  surface's 
microstructuie. 
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The  grass  surface  was  a piece  of  sod  with  fine,  densely  spaced  blades  of  grass.  The  blades  were  not 
small  enough  relative  to  the  laser-beam  diameter  to  avoid  a iarge  random  error  in  the  measurements.  Thus, 
only  the  statistical  properties  of  the  measurements  are  significant.  The  large  om  is  probably  due  mostly 

to  this  random  error.  However,  a significant  systematic  error  was  discernible.  Also,  from  the  absolute 
values  of  BR1DF  given  by  the  model,  the  directional-hemispherical  reflectance  at  normal  incidence  was 
estimated  to  be  of  the  order  of  60  percent,  which  is  inconsistent  with  the  S to  20  percent  measured  for 
several  samples  of  grassy  meadow.  However,  a value  of  e much  larger  than  unity  is  realistic;  that  is,  the 
surfaces  of  the  blades  of  grass  tend  to  be  vertical,  resulting  in  a concentration  of  the  microsurface  area 
at  high  angles  from  the  macrosurface,  and  the  model  irregularity  with  an  e greater  than  1 also  has  this 
property. 


For  all  the  surfaces,  the  optimization  gave  targe  values  for  k;  greater  than  0.3.  Such  values  occur 
only  for  metals  (reference  44).  and  these  surfaces  are  obviously  dielectric.  Even  a value  of  k of  only 

0.01  occurs  (reference  5)  only  for  an  extremely  highly  absorbing  dielectric  such  as  black  glass,  in  which 
95  percent  of  the  energy  is  absorbed  while  traveling  a distance  of  only  25  wavelengths.  In  one  paper 
(reference  48),  optimization  involving  k resulted  in  zero  values  for  dielectrics.  In  the  present  optimiza- 
tion (finding  the  values  of  the  parameters  giving  the  minimum  overall  deviation  between  theoretical  and 
experimental  data),  a trough  of  minima  was  found  in  the  space  of  the  parameters  n and  k.  For  each 
surface,  the  bottom  of  this  trough  was  nearly  level,  making  the  n k optimization  somewhat  uncertain. 
Also,  this  trough  intersected  the  k = 0 axis.  It  is  possible  that  the  position  of  the  minimum  found  along 
the  trough  for  each  surface  was  a random  occurrence  resulting  from  random  errors  in  the  experimental 
data  and  that  the  actual  minimum  was  at  k = 0.  For  one  surface,  the  3M  Black  Velvet,  two  minima 
were  found,  one  at  k = 0. 


Summarizing  the  comparison  of  the  model  to  the  measurements,  it  has  been  shown  that  the  com- 
parison is  in  general  reasonably  good  and  that  discrepancies  can  be  mostly  explained  by  the  existence  of 
gross  deviations  of  some  of  the  measured  surfaces  from  the  assumed  surface,  such  as  the  existence  of 
two  reflecting  surface  materials  instead  of  just  one  and  the  existence  of  a significant  directional  depend- 
ence of  the  surface  structure.  Since  most  of  the  data  were  taken  outside  the  plane  of  incidence,  these 
results  lend  to  verify,  outside  the  plane  of  incidence,  the  ability  of  interlace  plus  Lambertian  reflection 
to  give  most  of  the  reflection  m most  directions  for  most  surfaces.  Also,  it  has  been  shown  that  the 
ellipsoid  of  revolution  average  surface  irregularity  is  a useful  concept,  since  (1)  for  many  surfaces,  one 
can  visualize  the  ellipticily  of  the  representative  ellipsoid,  and  (2)  the  ellipsoid  gives  an  Improved  surface 
structure  I'm., Hon.  Last,  there  has  been  some  verification  of  the  correctness  of  the  normalization,  since 
the  optimization  was  done  on  relative  data,  and  the  results  usually  gave  absolute  directional -hemispher- 
ical reflectances  comparable  to  typical  measured  values. 


CONCLUSIONS 

The  following  significant  conclusions  may  Ik  drawn  about  the  ray-optics  theory  of  light  reflection 
from  a rough,  ait/matcnal  interface. 

1.  The  facet  representation  of  ihc  tough  interface  must  give  the  same  reflection  as  the  actual  inter- 
face. winch  is  made  up  of  curved  surfaces. 

2.  For  any  given  rough  surface,  there  exists  a single  optically  smooth  cutved  surface  of  revolution 
(average  irregularity)  of  very  restricted  shape  that  wilt  reflect  light  in  the  same  distribution  as 
that  reflected  by  the  rough  interface. 
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3.  Modeling  this  average  irregularity  as  an  ellipsoid  of  revolution  gives  a surface  structure  function 
that  is  much  more  accurate  and  useful  than  previously  existing  ones. 

4.  There  now  exists  a reflectance  model  that  can  be  normalized  giving  a reflectance-distribution 
function  that  is  absolute,  and  the  normalization  has  been  verified  experimentally. 

5.  The  combination  of  ray  reflection  from  the  rough  interface  plus  some  Lambertian  reflection 
applies  rather  well  outside  the  plane  of  incidence  as  well  as  within  it  for  a variety  of  commonly 
occurring  rough  surfaces. 
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APPENDIX 

PROOF  THAT  x IS  FINITE  WHEN 
THE  SLOPE  OF  h<x)  IS  VERTICAL 


The  contention  that  x is  still  finite  when  the  slope  of  h(x)  is  vertical  can  be  more  rigorously 
verified  by  the  following  derivation  of  the  upper  bound  x0  on  the  value  of  x at  which  h'(x)  = 
Equation  (7)  can  be  rewritten  as 


h"(x)  = -Cx  (I  +(hV|J'-/D(a>. 


which  upon  integration  gives  a ntonotonically  decreasing  function  h’(x)  provided  that  0<D(aK‘x>.  A ! 

lower  bound  on  l»"(x)l  is  given  from  i 

! 

! 


lt’’(x)l>Cxli,lJ/DV 

where  D*  is  the  upper  bound  on  IXa).  Since  both  h‘(x)  and  li“(xl  are  negative,  we  have 

h‘7(hV>Cx/l>*. 

which  upon  integration  give' 


l/lhjl“  < l/(hj  r ♦ C{.\‘  ♦ 

where  h’,  and  h,  ate  values  of  the  slope  h’lx)  at  v,  a.sd  a,,  respectively.  Finally.  It’,  E o* curs  at 
some  value  of  x , . where 


x,  = xn  < (xj  ♦ D*'C(h',)2|" 


litis  proves  the  contention,  since  the  expression  on  the  right  is  a finite  upper  bound  on  xfl. 
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